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Relativistic cosmological perturbation analyses can be made based on several different fundamen- 
tal gauge conditions. In the pressureless limit the variables in certain gauge conditions show the 
correct Newtonian behaviors. Considering the general curvature (K) and the cosmological constant 
(A) in the background medium, the perturbed density in the comoving gauge, and the perturbed 
velocity and the perturbed potential in the zero-shear gauge show the same behavior as the New- 
tonian ones in general scales. In the first part, we elaborate these Newtonian correspondences. In 
ON ' the second part, using the identified gauge- invariant variables with correct Newtonian correspon- 

ON . dences, we present the relativistic results with general pressures in the background and perturbation. 

We present the general super-sound-horizon scale solutions of the above mentioned variables valid 
for general K, A, and generally evolving equation of state. We show that, for vanishing K, the 
super-sound-horizon scale evolution is characterised by a conserved variable which is the perturbed 
three-space curvature in the comoving gauge. We also present equations for the multi-component 
hydrodynamic situation and for the rotation and gravitational wave. 
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NO . I. INTRODUCTION 

o 
o 

The analysis of gravitational instability in the expanding universe model was first presented by Lifshitz in 1946 
^\ , in a general relativistic context Q. Historically, the much simpler, and in hindsight, more intuitive Newtonian 
study followed later H. The pioneering study by Lifshitz is based on a gauge choice which is commonly called the 
synchronous gauge. As the later studies have shown, the synchronous gauge is only one way of fixing the gauge 
~ Q ■ freedom out of several available fundamental gauge conditions |5p||. As will be summarized in the following, out of 
the several gauge conditions only the synchronous gauge fails to fix the gauge mode completely and thus often needs 
more involved algebra. As long as one is careful of the algebra this gauge choice does not cause any kind of intrinsic 
problem; there exist, however, some persisting algebraic errors in the literature, see Sec. The gauge condition 
which turns out to be especially suitable for handling the perturbed density is the comoving gauge used first by Nariai 
in 1969 Q. Since the comoving gauge condition completely fixes the gauge transformation property, the variables in 
this gauge can be equivalently considered as gauge-invariant ones. As mentioned, there exist several such fundamental 
gauge conditions each of which completely fixes the gauge transformation properties. One of such gauge conditions 
suitable for handling the gravitational potential and the velocity perturbations is the zero-shear gauge used first by 
Harrison in 1967 The variables in such gauge conditions are equivalently gauge-invariant. Using the gauge freedom 
as an advantage for handling the problem was emphasized by Bardeen in 1988 ||]. In order to use the ga uge c hoice 
as an advantage a gauge ready method was proposed in || which will be adopted in the following, see Sec. [[ib|. 

The variables which characterize the self gravitating Newtonian fluid flow are the density, the velocity and the 
gravitational potential (the pressure is specified by an equation of state), whereas, the relativistic flow may be 
characterized by various components of the metric (and consequent curvatures) and the energy- momentum tensor. 
Since the relativistic gravity theory is a constrained system we have the freedom of imposing certain conditions on 
the metric or the energy-momentum tensor as coordinate conditions. In the perturbation analyses the freedom arises 
because we need to introduce a fictitious background system in order to describe the physical perturbed system. The 
correspondence of a given spacetime point between the perturbed spacetime and the ficticious background one could 
be made with certain degrees of freedom. This freedom can be fixed by suitable conditions (the gauge conditions) 
based on the spacetime coordinate transformation. Studies in the literature show that a certain variable in a certain 
gauge condition correctly reproduces the corresponding Newtonian behavior. Although the perturbed density in 
the comoving gauge shows the Newtonian behavior, the perturbed potential and the perturbed velocity in the same 
gauge do not behave like the Newtonian ones; for example, in the comoving gauge the perturbed velocity vanishes by 
the coordinate (gauge) condition. It is known that both the perturbed potential and the perturbed velocity in the 
zero-shear gauge correctly behave like the corresponding Newtonian variables || . 
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In the first part (Sec. til) we will elaborate establishing such correspondences between the relativistic and Newtonian 
perturbed variables. Our previous work on this subject is presented in JT^JTTJ] ; compared with |l(J in this work we 
will explicitly compare the relativistic equations for the perturbed density, potential and velocity variables in several 
available gauge conditions with the ones in the Newtonian system. We will include both the background spatial 
curvature (K) and the cosmological constant (A). In the second part (Sec. |Tv| ), using the variables with correct 
Newtonian correspondences, we will extend our relativistic results to the situations with general pressures in the 
background and perturbations. We will present the relativistic equations satisfied by the gauge-invariant variables 
and will derive the general solutions valid in the super-sound- horizon scale (i.e., larger than Jeans scale) considering 
both K and A, and the generally evolving ideal fluid equation of state p = p(ji). 

Section H is a review of our previous work displaying the basic equations in both Newtonian and relativistic contexts 
and summarizing our strategy for handling the equations. In Sec. [II we consider a pressureless limit of the relativistic 
equations. We derive the equations for the perturbed density, the perturbed potential and the perturbed velocity in 
several different fundamental gauge conditions. By comparing these relativistic equations in several gauges with the 



Newtonian ones we identi 
in general scales. In Sec. 



the gauge conditions which reproduce the correct Newtonian behavior for certain variables 
] we present the fully relativistic equations for the identified gauge-invariant variables with 
correct Newtonian limits, now, considering the general pressures in the background and perturbations. We derive 
the general large-scale solutions valid for general K and A in an ideal fluid medium, which are thus valid for general 
equation of state of the form p = p(fj.), but with negligible entropic and anisotropic pressures. The solutions are 
valid in the super-sound-horizon scale, thus are valid effectively in all scales in the matter dominated era where the 
sound-horizon (Jeans scale) is negligible. We discuss several quantities which are conserved in the large-scale under the 
changing background equation of state. These are the perturbed three-space curvature in several gauge conditions, and 
in particular, we find the three-space curvature perturbation in the comoving gauge shows a distinguished behavior: 
for K = (but for general A), it is conserved in a super-sound- horizon scale independently of the changing equation 
of state. For completeness we also summari ze th e case with multiple fluid components in Sec 
rotation and the gravitational wave in Sec. [VD. Sec. IV is the highlight of the present work. 
We set c = 1. 
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II. BASIC EQUATIONS AND STRATEGY 



A. Newtonian Cosmological Perturbations 



Since the Newtonian perturbation analysis in the expanding medium is well known we begin by summarizing the 
basic equations Ilfljl2,0 . The background evolution is governed by 
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H = -^-g 2 + — , goto, 
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(1) 



where we allowed the general curvature (total energy) and the cosmological constant; a(t) is a cosmic scale factor, 
H(t) = a/a, and g(t) is the mass density. In Newtonian theory the cosmological constant can be introduced in the 
Poisson equation by hand as V 2 $ = 4irGg — A. Perturbed parts of the mass conservation, the momentum conservation 
and the Poisson's equations are [see Eqs. (43,46) in [Eof]: 



Sg + 3H5g = gSv, 

a 



k k 2 

Sv + HSv = -(5$. -5$ = 4nGSg, 

a a z 



(2) 



where Sg(k, t), Sv(k,t) and S$(k, i) are the Fourier modes of the perturbed mass density, velocity, and gravitational 
potential, respectively [see Eqs. (44,45) in |10|]; k is a comoving wave number with V 2 — > ~k 2 /a 2 . [For linear 
perturbations the same forms of equations are valid in the configuration and the Fourier spaces. Thus, without 
causing any confusion, we often ignore distinguishing the Fourier space from the configuration space by an additional 
subindex.] Equation (g) can be arranged into the closed form equations for 5 (= Sg/g), Sv and <5<I> as: 



S + 2HS - inGgS 
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Sv + 3HSv + { H + 2H 2 - 4nGg) Sv = 0, 
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We note that Eqs. are valid for general K and A. Although redundant, we presented these equations for later 

comparison with the relativistic results. The general solutions for S, Sv and 5<S> immediately follow as (see also Table 

i in my. 



<f(k,t) = k 2 
5v(k,t) = - 

<J$(M) = - 
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The coefficients C(k) and d(k) are two integration constants, and indicate the relatively growing and decaying so- 
lutions, respectivel y; the coefficients are matched in accordance with the solutions with a general pressure in Eqs. 
(|l7[jls| ) using Eq. ( p9| ) (different dummy variables for t in the integrands are assumed). In the literature we can find 
various analytic form expressions of the above solutions in diverse situations, see . 



B. Relativistic Cosmological Perturbations 

In the relativistic theory the fundamental variables are the metric and the energy-momentum tensor. As a way of 
summarizing our notation we present our convention of the metric and the energy-momentum tensor. As the metric 
we consider a spatially homogeneous and isotropic spacetime with the most general perturbation 

ds 2 = -a 2 (1 + 2a) drj 2 - a 2 + B a ) d V dx a + a 2 [ 5 g (1 + 2<p) + 2 7 | Q/3 + 2C W/3) + 2C Q/3 ] dx a dxf>, (9) 

(3) 

where = r\ with dt = adrj, Indices a, p . . . are spatial ones and g a i is the three-space metric of the homogeneous and 
isotropic space. a(x,t), (3(x,t), j(x,t), and <p(x,t) indicate the scalar-type structure with four degrees of freedom. 
The transverse B a (x,t) and C a (x,t) indicate the vector- type structure with four degrees of freedom. The transverse- 
tracefree C a p (x, t) indicates the tensor-type structure with two degrees of freedom. The energy-momentum tensor is 
decomposed as: 



T ° = -ft - Sfi, T° = -j(ll + p) v a + 8TW, 
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v i v /3 -^A^j ( T + ( 5r (t ' ) ; + ( 5r ( * ) ;, (io) 

where /i(k, t) = p,(t) + 5//(k, t) and p = p + Sp are the energy density and the pressure, respectively; an overbar 
indicates a background order quantity and will be ignored unless necessary, v is a frame-independent velocity related 
perturbed order variable and a is an anisotropic pressure. Situation with multiple fluids will be considered in Sec. 



IV C. 5T a, ST g, and 5T ^ are vector and tensor-type perturbed order energy-momentum tensor. All the 



(3) 

spatial tensor indices and the vertical bar are based on g y a + as the metric. The three types of structures are related 
to the density condensation, the rotation, and the gravitational wave, respectively. Due to the high symmetry in the 
background the three types of perturbations decouple from each other to the linear order. Since these three types 
of structures evolve independe ntly, we can handle them separately. Evolutions of the rotation and the gravitational 



wave will be discussed in Sec. IV D 



For the background equations we have [Eq. (6) in |lq]]: 

H 2 = ^ fl _K + A h = -AkG{»+p) + - 2 , A=-3^(M + P)- (11) 

These equations follow from Gq and G„ — 3Gq components of Einstein equations and Tq.^ = 0, respectively; the third 
equation follows from the first two equations. For p = and replacing p with g Eq. (|ll|) reduces to Eq. ([!]). The 
relativistic cosmological perturbation equations without fixing the temporal gauge condition, thus in a gauge ready 
form, were derived in Eq. (22-28) of M [see also Eqs. (8-14) in |TJ]; we use 5 = Sfi/fi and v = — (k /(£)"$>/ ((j, + p)]: 
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3 {-y + Ha) + k —x, (12) 



a 2 ' 
k 2 — "\K 

5 tp + Hk = -AirGfiS, (13) 

gt 

k 2 — 3K a 
k — X = 12wG(fi + p)-v, (14) 

X + H X -a-<p = 8ttGct, (15) 
k + 2Hk = ( -J - 3i?) a + 4ttG (l + 3c 2 ) fj.5 + UirGe, (16) 



5 + 'SH (c 2 - w) 5 + 3H~ = (1 + w) f K - 3Ha - '^-vj , (17) 

/ 2 \ k k ( 9r e 2 k 2 - 3K a\ 

V+(l- 3c 2 s ) Hv = -a+ — [c 2 s 5+-- - 5 , 18 

v sj a a (1 + w) \ (i 3 a 2 (i J v ' 

where we introduced a spatially gauge-invariant combination \ = a(/3 + a-y). According to their origins, Eqs. ( p"2| - 
[l8| ) can be called: the definition of k, ADM energy constraint (Gg component of Einstein equations), momentum 
constraint (G" component), ADM propagation [G% — ^SgG^ component), Raychaudhuri (G" — Gq component), 
energy conservation (T^. b = 0), and momentum conservation (T^. b = 0) equations, respectively. The perturbed metric 
variables <p(k, t), k(\c, t), x(k, t) and a(k, t) are the perturbed part of the three space curvature, expansion, shear, and 
the lapse function, respectively; meanings for ip, k, and x are based on the normal frame vector field. The isotropic 
pressure is decomposed as 

5p(k,t)=c 2 s (t)5[i(k,t)+e{k,t), c 2 s = ^-, w(t) = -. (19) 

H n 

Under the gauge transformation x a = x a + £ a , the variables transform as (see Sec. 2.2 in ||): 
a = a-^, <p = <p-H?, x = X-C\ « = « + ( 3H - Kr) ? , 
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k 



V = V 



a 



f, 5 = 6 + 3(1 +w) H?, (20) 



where £ = e 8 ^, and e and cr are gauge- invariant. [Due to the spatial homogeneity of the background, the effect from 
the spatial gauge transformation has been trivially handled; x an d v are the spatially gauge-invariant combinations of 
the variables, and the other metric and fluid variables are naturally spatially gauge-invariant 0.] Thus, the perturbed 
variables in Eqs. (|l2|-p"8"|) are designed so that any one of the following conditions can be used to fix the freedom based 
on the temporal gauge transformation: a = (the synchronous gauge), if = (the uniform-curvature gauge), x = 
(the zero-shear gauge), k = (the uniform-expansion gauge), v/k = (the comoving gauge), and 5 = (the uniform- 
density gauge). These gauge conditions contain most of the gauge conditions used in the literature concerning the 
cosmological hydrodynamic perturbation in a single component situation; for the multiple component case see Eq. 
(|69|). Of course, we can make an infinite number of different linear combinations of these gauge conditions each of 
which is also suitable for the temporal gauge fixing condition. Our reason for choosing these conditions as fundamental 
is partly based on conventional use, but apparently, as the names of the conditions indicate, the gauge conditions also 
have reasonable meanings based on the metric and the energy-momentum tensors. 

Each one of these gauge conditions, except for the synchronous gauge, fixes the temporal gauge transformation 
property completely. Thus, each variable in these five gauge conditions is equivalent to a corresponding gauge- 
invariant combination (of the variable concerned and the variable used in the gauge condition). We proposed a 
convenient way of writing the gauge-invariant variables For example, we let 

aH aH k k 
5 V = 5+ 3(1 +w)— v = 3(1 + w)— vg, tp x = f - Hx = -Hx v , v x =v x = Xv (21) 

The variables 5 V , v x and ip x are gauge- invariant combinations; 5 V becomes 5 in the comoving gauge which takes 
v/k = as the gauge condition, etc. In this manner we can systematically construct the corresponding gauge- 
invariant combination for any variable based on a gauge condition which fixes the temporal gauge transformation 
property completely. One variable evaluated in different gauges can be considered as different variables, and they 
show different behaviors in general. In this sense, except for the synchronous gauge, the variables in the rest of five 
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gauges can be considered as the gauge-invariant variables. [The variables with a subindex a are not gauge-invariant, 
because those are equivalent to variables in the synchronous gauge.] Although 5 V is a gauge-invariant variable which 
becomes S in the comoving gauge, S v itself can be evaluated in any gauge with the same value. 

Equations ( |l2[jl^ ) are the basic set of gauge-ready form equations for hydrodynamic perturbations proposed in 
M: "The moral is that one should work in the gauge that is mathematically most convenient for the problem at 
hand", §. Accordingly, in order to use the available gauge conditions as the advantage, Eqs. (J12 - 18 ) are designed 
for easy implementation of the fundamental gauge conditions. Using this gauge-ready formulation, complete sets of 
solutions for the six different gauge conditions are presented in an ideal fluid case 0], and in a pressureless medium 
p0| . Equations ( pT[]l9| ) include general pressures which may account for the nonequilibrium or dissipative effects 
in hydrodynamic flows in the cosmological context with general K and A. The equations are expressed in general 
forms so that the fluid quantities can represent successfully the effects of the scalar field and classes of generalized 
gravity theories jg]. The applications of the gauge-ready formalism to the minimally coupled scalar field and to the 
generalized gravity theories have been made in p6[ . 



III. NEWTONIAN CORRESPONDENCE OF THE RELATIVISTIC ANALYSES 



In 



[10|,|11| we developed arguments on the correspondences between the Newtonian and the relativistic analyses. 
In order to reinforce the Newtonian correspondences of certain (gauge-invariant) variables in certain gauges, in the 
following we will present the closed form differential equations for <5, v and (p in the six different gauge conditions and 
compare them with the Newtonian ones in Eqs. 

There are reasons why we should know about possible Newtonian behaviors of the relativistic variables. Variables 
in the relativitistic gravity are only parameters appearing in the metric and energy-momentum tensor. For example, 
later we will find that only in the zero-shear gauge does the perturbed curvature variable behave as the perturbed 
Newtonian potential, which we have some experience of in Newtonian physics. The same variable evaluated in other 
gauge conditions is simply other variable and we cannot regard it as the perturbed potential. As we have introduced 
several different fundamental gauge conditions it became necessary to identify the correct gauge where a variable shows 
the corresponding Newtonian behavior. As will be shown in this section, we rarely find Newtonian correspondences. In 
fact, there exists an almost unique gauge conditio n for e ach relativistic perturbation variable which shows Newtonian 



behavior. The results will be summarized in Sec. Ill G 



In this section we consider a pressureless medium, p — (thus w = = c 2 s ) and e = = a. However, we consider 
general K and A. 



A. Comoving Gauge 

As the gauge condition we set v/k = 0. Equivalently, we can set v/k = and leave the other variables as the 
gauge-invariant combinations with subindices v/k or simply v. From Eq. ( |l8| ) we have a v — 0. Thus the comoving 
gauge is a case of the synchronous gauge; this is true only in a pressureless situation. From Eqs. (p~^-|l8|) we can 
derive: 

S v + 2H5 V - 4irGfiS v = 0, (22) 
K 

Cp v + 3H(p v *<Pv = 0. (23) 

cr 

Thus, Eq. (^2|) has a form identical to Eq. (||). However, Eq. (|23|) differs from Eq. (||), and apparently we do not 
have an equation for v. 

For K = Eq. ( p3| ) leads to two solutions which are ip v oc constant and J* a~ 3 dt. Since we have a v = 0, for K = 0, 
from Eqs. ( |T^Jl4| ) we have tp v = 0. This implies that, for K = 0, the second solution, ip v oc J* a~ 3 dt, should have the 



vanishing coefficient, see Eq. (56). The solution of Eq. (|23|) for general K will be presented later in Eq. (53) 



B. Synchronous Gauge 

We let a = 0. This gauge condition does not fix the temporal gauge transformation property completely. Thus, 
although we can still indicate the variables in this gauge condition using subindices a without ambiguity, these 
variables are not gauge-invariant; see the discussion after Eq. (|52]). Equation (|l8|) leads to 
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Cg — , 



(24) 



which is a pure gauge mode. Thus, fixing c g = exactly corresponds to the comoving gauge. We can show that the 
following two equations are not affected by the remaining gauge mode. 



S a + 2HS a - 4irGn5 a = 0, 
K 

tp a + 3Hip a -ip a = 0. 



(25) 
(26) 



Equation (|25|) is identical to Eq. (|3j). This is because in a pressureless medium the behavior of the gauge modes 
happen to coincide with the behavior of the decaying physical solutions for 5 a and tp a . However, for the variables K a 
and \a the gauge mode contribution appears explicitly. 

Thus, in a pressureless medium, variables in the synchronous gauge behave the same as the ones in the comoving 
gauge, except for the additional gauge modes which appear in the synchronous gauge for some variables. In a 
pressureless medium, we can simultaneously impose both the comoving gauge and the synchronous gauge conditions. 
However, this is possible only in a pressureless medium, see the Appendix. In Sec. [v| we indicate some common errors 
in the literature based on the synchronous gauge. 



C. Zero-shear Gauge 

We let x = 0, and substitute the other variables into the gauge-invariant combinations with subindices X- We can 
derive: 



2k 2 /a 2 - 36ttG/x 1 + (H 2 - 2H)a 2 / (k 2 ~ 3K) 



-4irGn 
3Hv 



k 2 /a 2 - UirGn [1 + 3H 2 a 2 / (k 2 - 3K)] 

2 



H [L - 12-kGh 



1 - 3(H 2 - 2H)— 

1 'k 2 -3K 



5 X = 0, 



v x -r on v x -r[H+ 2H 2 - 4itGh) v x = 0, 
(p x + 4Hlp x + (h+ 3H 2 - 4nGv) ip x = 0. 



(27) 

(28) 
(29) 



Thus, Eqs. ( p8| , p9| ) have forms identical to Eqs. (||,|]), respectively. However, only in the small-scale limit is the 
behavior of S x the same as the Newtonian one. 

D. Uniform-expansion Gauge 

We let k = 0, and substitute the remaining variables into the gauge-invariant combinations with subindices k. We 
can derive: 



S K 



2H - _ ^ . ] S K 



2H 



k 2 /a 2 + 6H 2 



4:TiGn 
AnG^i I k 2 /a 



k 2 /a 2 - 3H 

UttG^H 
k 2 /a 2 -oH ^ k 2 /a 2 -3H 

k 2 /a 2 + 3H 2 



k 2 /a 2 - 3H 



2 

2 /„2 



3H 



H + H z - AttG^i 



k 2 /a 2 - 3H 



H 



4ttGh 



\2ttGhH 

k 2 /a 2 - 3H 



k 2 /a 2 - 3H 



0. 



(30) 



AH 



UirGfiH 

k 2 /a 2 - 3H 



k 2 /a 2 -3H \ 47rG*Ai 
k 2 /a 2 + 9H 2 



H + 3H Z - 4:TrGn- 



k 2 /a 2 ~3H \k 2 /a 2 -3H 



v K = 0, 
12-kGijlH 



0. 



(31) 
(32) 



In the small-scale limit we can show that Eqs. (|3JJ-|32|) reduce to Eqs. (|3jj5J), respectively. Thus, in the small-scale 
limit, all three variables S K , v K and (p K correctly reproduce the Newtonian behavior. However, outside or near the 
(visual) horizon scale, the behaviors of all these variables strongly deviate from the Newtonian ones. 



G 



In Sec. 84 of |13[ we find that in order to get the usual Newtonian equations a coordinate transformation is made 
so that we have h = in the new coordinate. We can show that h = 2k in our notation. Thus, the new coordinate 
used in |l3| is the uniform-expansion gaugef]. 

E. Uniform-curvature Gauge 

We let ip = 0, and substitute the other variables into the gauge-invariant combinations with subindices (p. We have 



■■ {k 2 -3K)/a 2 + l^G l i - AnG^/a 2 , _ n 

v + (fc 2 - 3A)/a 2 + 12ttG/ v (fc 2 - 3A)/a 2 + UirGfi v ~ ' 



»», + + v v + (zH + 4i? 2 ) v v = 0. (34) 

In the small-scale limit Eq. ( |33] ) reduces to Eq. (^). In the uniform-curvature gauge, the perturbed potential is set to 
be equal to zero by the gauge condition. The uniform-curvature gauge condition is known to have distinct properties 
in handling the scalar field or the dilaton field which appears in a broad class of the generalized gravity theories 

MM- 

F. Uniform-density Gauge 

We let 5 = 0, and substitute the other variables into the gauge-invariant combinations with subindices S. We have 

v 5 + 2 (2H + ^ J v 5 + 3 (h 2 + Hj v s = 0, (35) 
.. (P-SKy^ + mnG^ . A^k 2 /a 2 

Vs + 2 V ~ 3X)/a 2 + l 27 rG^ " (fc 2 - 3A)/a 2 + 12nG^ S = °" (36) 

These equations differ from Eqs. (|],|]). Of course, we have no equation for 5 which is set to be equal to zero by our 
choice of the gauge condition. 

G. Newtonian Correspondences 

After a thorough comparison made in this section we found that equations for S in the comoving gauge (S v ), and 
for v and tp in the zero-shear gauge (v x and ip x ) show the same forms as the corresponding Newtonian equations in 
general scales and considering general K and A. Using these gauge-invariant combinations in Eq. (pl|), Eqs. ( p"2]ji"i^ ) 
can be combined to give: 

• k 2 -3K k k 2 -3K 

+ Hv x = ip x , ip x = A-kG\io v . (37) 



7 ' V X' u 

aK (x lij 

Comparing Eq. ( ^7| ) with Eq. (Jfy we can identify either one of the following correspondences: 



k 2 -3K „ k 2 -3K 

' I Newtonian •> 

k 2 



S v ^Newtonian, p V x <-> Sv | Newtonian , p fx ^ | Newtonian, (38) 



k 2 - 3K 



S v <-> (^Newtonian, V X ^ | Newtonian , — f X I Newtonian- (39) 



x This was incorrectly pointed out after Eq. (49) in |lo| : In JToJ it was wishfully mentioned that in Sec. 84 of Jl^] the gauge 
transformations were made into the comoving gauge for 8 and into the zero-shear gauge for v and <p, respectively. 
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In [|10| we proposed the correspondences in Eq. (|3£|), but the ones in Eq. ( p9[ ) also work well. In fact, the gravitational 
potential identified in Eq. (|3^) is the one often found in the Newtonian limit of the general relativity, e.g. in [[l2| . 
Using Eqs. 



mi 



37) we can also identify the following relations: 



S v = 3H-vs, 
k 



uk 



-Xv 



k 2 - 3K 



fx 



-HXip- 



(40) 



From a given solution known in one gauge we can derive all the rest of the solutions even in other gauge conditions. 
This can be done either by using the gauge-invariant combination of variables or directly through gauge transforma- 
tions. The complete set of exact solutions in a pressureless medium is presented in | jlC| |. From our study in this section 
and using the complete solutions presented in Tables 1 and 2 of |hJ we can identify variables in certain gauges which 
correspond to the Newtonian ones. These are summarized in Table 1. 

[[TABLE 1.]] 

As mentioned earlier, all three variables in the uniform-expansion gauge show Newtonian correspondence in the small- 
scale; however all of them change the behaviors near and above the horizon scale. In the small-scale limit, except for 
the uniform-density gauge where 5 = 0, S in all gauge conditions behaves in the same way ||. Meanwhile, since only 
5 V (6 in the comoving gauge) shows the Newtonian behavior in a general scale, it may be natural to identify S v as the 
one most closely resembling the Newtonian one. 

Notice that, although we have horizons in the relativistic analysis the equations for S v , v x and <p x keep the same 
form as the corresponding Newtonian equations in the general scale. Considering this as an additional point we regard 
S Vl v x and <p x as the one most closely corresponding to the Newtonian variables. This argument will become stronger 
as we consider the case with a general pressure in the next section. 



IV. RELATIVISTIC COSMOLOGICAL HYDRODYNAMICS 



A. General Equations and Large-scale Solutions 



In the previous sections we have shown that the gauge- invariant combinations S v , v x and ip x behave most similarly 
to the Newtonian S = Sg/ g, 5v and 5<J>, respectively. The equations remain the same in a general scale which includes 
the superhorizon scales in the relativistic situation considering general K and A. In this section, we will present the 
fully general relativistic equations for S v , v x and (p x including the effects of the general pressure terms. Equations 
( |l2[ - |lq ) are the basic set of perturbation equations in a gauge-ready form. 

From Eqs. Eqs. (f|JI|), Eqs. (|l|,|l|), and Eqs. (0,0,[l|) , we have, respectively: 



k k 2 

S v - 3HwS v = - (1 + w) 

a k z 



v x + Hv x = — LD X H r 

x x a x a(l+w) 



3A 



2H- 



3Ka 



c 2 s S v + -- 8ttG (1 + w) a 



2 k 2 - 3K a 

3 a 2 /i 



ip x = AirGfidy, 

n z 



fx 



Hip x = -4vrG (n + p) -v x - 8nGHa. 



(41) 

(42) 

(43) 
(44) 



Considering either one of the correspondences in Eqs. (^8|,^9|) we immediately see that Eqs. ( fil"|{4~4| ) have the correct 
Newtonian limit expressed in Eqs. (|,||). Equations (glfg) were presented in H . 

Combining Eqs. (j4l|-|44"|) we can derive closed form expressions for the (5„ and ip x which are the relativistic counterpart 
of Eqs. (|,|). We have: 



3c 2 s 
1 



6w) HS V 4 
H 2 



a 2 H 
k 2 - 



K 



4nGn (1 - 6c 2 s + 8w - 3w 2 ) + 12 (w - c 2 s ) — + (3c 2 s - 5w) A 



a 3 A* , 



a(ji +p) V H 

+■ 2Ha + 2 



-3K 1 
— < ' 



+ c-S v 

-\ k ^ + lH + 3(l + c 2 )H 2 



(45) 
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Cp x + (4 + 3c 2 ) H<p x + 

h+p r # 2 



^ + 8?rG ( c 2 _ w ) _ 2 (1 + 3c 2 ) 4 



A 



<Px 



H 



a(fj,+p) 



stresses. 



(46) 



These two equations are related by Eq. (|43]_which resembles the Poisson equation. Notice the remarkably compact 
forms presented in the second steps of Eqs. (4^j46]). It may be an interesting exercise to show that the above equations 
are indeed valid for general K and A, and for the general equation of state p = p{p); use w = — 3H(1 + w)(Cg — w). 
Equation (|45l ) became widely known through Bardeen's seminal paper in B . In a less general context but originally 
Eqs. ( [45| , [46D were derived by Nariai B and Harrison ||, respectively; however, the compact expressions are new 
results in this paper. 



If we ignore the entropic and anisotropic pressures (e = = a) on scales larger than Jeans scale Eq. (46) immediately 
leads to a general integral form solution as 



Vx (k,t)=4nGC(k)- f 
a Jo 

= C(k) 



H 



a(y + p) 
H 2 

1 - 



H 



dt 



d(k) 



H 



d(k). 



(47) 



This solution was first derived in Eq. (108) of Jl7]]0, see also Eq. (55) in 
(E3LE4I), respectively, as 



Solutions for 6 V and v x follow from Eqs. 



5 v (k,t) 
v x (k,t) 



k 2 



— 2 ^x( k > 0' 



4ttG(^ + p)a 2 



C(k) 



K 
a 



H I a 

'0 



K 

A2 



dt 



Hd(k) 



(48) 
(49) 



We stress that these large-scale asymptotic solutions are valid for the general K , A, and p — p(p); C(k) and d(k) are 
integration constants considering the general evolution of p — p(p). We also emphasize that the large-scale criterion 
is the sound-horizon, and thus in the matter dominated era the solutions are valid even far inside the (visual) horizon 

correctly 



as long as the scales are larger than Jeans scale. In the pressureless limit, considering Eq. (p9[), Eqs. ( [47 - 
reproduce Eqs. (^-g)- 

Using the set of gauge-ready form equations in Eqs. (|l^jl^ ) the complete set of corresponding general solutions for 
all variables in all different gauges can be easily derived; for such sets of solutions with less general assumptions see 
@|5|. For K = = A and w = constant we have a cx t 2 /[ 3 ( 1+w )l and Eqs. @^|) become: 



MM) 



3(1 



5 + 3w 



-C(k) 



1 



3(1 



d(k) 



oc constant, 

2(l + 3w) 

8 v (k,t) oc 1 3 < 1 +'») , 



v x (k, t) tx t~ 



5 + 3w 
t 3(l + w) 



w) at 

constant, 



a 



l+3w 



a 



a 

2 



a 



These solutions were presented in 
well known Newtonian behaviors. 



and in less general contexts but originally in pH 



(50) 

For w = we recover the 



B. Conserved Quantities 

In an ideal fluid, the curvature fluctuations in several gauge conditions are known to be conserved in the superhorizon 
scale independently of the changes in the background equation of state. From Eqs. (41,73) in jl5| we find [for 
K = = A, but for general p(n)]- 



2 We correct a typographical error in Eq. (129) of j^: co/a should be replaced by cq. 
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<Pv = ¥>6 = <Pk = <f a = C(k), (51) 
and the dominating decaying solutions vanish (or, are cancelled). The combinations follow from Eq. (po|) as: 



aH 



tp v = tp - 



-v, ips = (f- 



3(1 +w) 



H 



3H-k 2 /d- 



-k, tp a — tp — H adt. 



(52) 



The tp a combination which is tp in the synchronous gauge is not gauge- invariant; the lower bound of integration gives 
the behavior of the gauge mode; in Eq. ( |5l"| ) we ignored the gauge mode. The large-scale conserved quantity, £, 
proposed in [^,^8) is tps. In jig ] the above conservation properties are shown assuming K = = A; in such a case 
the integral form solutions for a complete set of variables are presented in Tabic 8 in p5| . In Eqs. ( ^7|^9| ) we have 
the large-scale integral form solutions in the case of general K and A. Thus, now we can see the behavior of these 
variables considering the general K and A. 

Evaluating tp v in Eq. ( p2| ) in the zero-shear gauge, thus tp v = tp x — (aH/k)v x , and using the solutions in Eqs. 
([l7],[l9|) we can derive 



b 



<^(k,t) = C(k){l + 
= C(k) 1 



K 



1 



H 

1 

a 



a 2 Ai:G(fi+p) 
K H/a [* a{[i+p) 
H 2 



a M + P J a 

Thus, for K = (but for general A) we have 



dt 



K 



H/a 



K H/a 
a 2 4irG(fi+p) 



d(k) 



a 2 4ttG(^ + p) 



d(k). 



^(k,t) = C(k), 



(53) 



(54) 



with the vanishing decaying solution. The disappearance of the decaying solution in Eq. (|5J) implies that the 
dominating decaying solutions in Eqs. @||) cancel out for K = 0. In fact, for K = 0, from Eqs. (|^,[l|,[l^,|l^,[l8|) 

we can derive 



r a 3 (pt+p) . 
c 2 H 2 ^ 



-tp v = stresses. 



For a pressureless case (c 2 — 0), instead of Eq. ( |55| ) we have W = 0; see after Eq. 
ignoring the stresses, we have an integral form solution 



<p v = C(x) - ^(x) 



c 2 H 2 



o a 3 (fi + p) 



-dt. 



(55) 

In the large-scale limit, and 
(56) 



Compared with the solution in Eq. (54), the D term in Eq. ( J56| ) is c s (k/aH) order higher than d terms in the other 
gauge. Therefore, for K — 0, tp v is conserved for the generally varying background equation of state, i.e., for general 
p = p(jj). Since the solutions in Eqs. (^?]|l8|,[l9|) are valid in super-sound-horizon scale, the conservation property in 
Eq. ( |54j ) is valid in all scales in the matter dominated era. 
For tps, by evaluating it in the comoving gauge we have 



Using Eqs. (M,f48h we have 



tp 5 (k,t) =C(k) 



\ a 2 12nG(ii+p) 



tps =tp v + 



1-* 

a 



3(1 + w)' 



,/ j 1 - — ) «// 



H/a 



a 2 l2nG(fi + p) 



d(k). 



(57) 



(58) 



Since the higher order term ignored in Eq. ( |53| ) is c 2 (k/aH) 2 order higher, in the medium with c s ~ 1(= c), the terms 
involving (fc/a) 2 /(/i + p) are not necessarily valid. Thus, to the leading order in the superhorizon scale we have 



¥>i(k,t)=C(k). 



(59) 



Thus, in the superhorizon scale tps is conserved apparently considering K and A. [In the case of nonvanishing K we 
may need to be careful in defining the superhorizon scale. Here, as the superhorizon condition we simply took k 2 term 
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to be negligible. In a hyperbolic (negatively curved) space with K = — 1, since i?(3) = 6K/a 2 , the distance a/y/\K~\ 
introduces a curvature scale: on smaller scales the space is near flat and on larger scales the curvature term dominates. 
In the unit of the curvature scale k 2 less than one (~ \K\) corresponds to the super-curvature scale, and k 2 larger than 
one corresponds to the sub-curvature scale. Thus, when we ignored the k 2 term compared with \K\ in a negatively 
curved space, we are considering the super-curvature scale. In other words, for the hyperbolic backgrounds k 2 takes 
continuous value with k 2 > where k 2 > 1 corresponds to subcurvature mode whereas < k 2 < 1 corresponds to 
supercurvature mode, and our large-scale limit took k 2 — > 0. A useful study in the hyperbolic situation can be found 
in psfl -l However, in a medium with the negligible sound speed (like the matter dominated era), Eq. (f58| ) is valid 
considering the k 2 terms. Thus, as the scale enters the horizon in the matter dominated era <ps wildly changes its 
behavior and is no longer conserved inside the horizon; inside the horizon ip$ is dominated by the S part which shows 
Newtonian behavior of S in most of the gauge conditions, see Table 1. 

For ip K , by evaluating it in the uniform-density gauge and using Eq. (|l3|), which gives k$ = (k 2 — 3K)/(a 2 H)ips, 
we have 



3K 



l2nG{fi + p)a 2 



(60) 



Thus, in the superhorizon scale and for K = 0, tp K is conserved. As in tps, in a medium with the negligible sound 
speed, the solution in Eq. ( |60| ) with Eq. ( p8| ) is valid considering the k 2 terms. Thus, as the scale enters the horizon 
in the matter dominated era tp K also changes its behavior and is no longer conserved inside the horizon. 

Now we summarize: In the superhorizon scale tp v and tp K are conserved for K = 0, while tps is conserved considering 
the general K. In the matter dominated era, for K = 0, tp v is still conserved in the super-sound-horizon scale which 
virtually covers all scales, whereas tps and ip K change their behaviors near the horizon crossing epoch and are no longer 
conserved inside the horizon. In this regard, we may say, ip v is the best conserved quantity. Curiously, although tp x is 
the variable most closely resembling the Newtonian gravitational potential, it is not conserved for changing equation 
of state, see Eq. (fT?]); however, it is conserved independently of the horizon crossing in the matter dominated era for 
K = = A, see Eq. ©. 

Similar conservation properties of the curvature variable in certain gauge conditions remain true for models based 
on a minimally coupled scalar field and even on classes of generalized gravity theories. In the generalized gravity 
the uniform-field gauge is more suitable for handling the conservation property, and the uniform-field gauge coincides 
with the comoving gauge in the minimally coupled scalar field. Thorough studies of the minimally coupled scalar field 
and the generalized gravity, assuming K = = A, are made in p^pOj ] and summarized in pi] ]. 

C. Multi-component situation 

We consider the energy-momentum tensor composed of multiple components as 

Tab = ^ T {l)ab> T (i)a-b = Q(i)o) zl Q(l)a = °! ( 61 ) 
(0 (0 

where i, l(= 1,2, ... ,n) indicates n fluid components, and Q{i) a considers possible mutual interactions among fluids. 
Since we are considering the scalar-type perturbation we decompose the Q(i) a in the following way: 

Q(i)0 = ~ a (l + a )Q(i)i Q(i) = Q{%) + 8Q(i)> Q(i)a = J(i),a- (62) 



The scalar-type fluid quantities in Eq. (|10|) can be considered as the collective ones which are related to the individual 
component as: 

(0 (0 (0 (0 

(ji + p)v = 5^(M(0 +P(i))v(i), cr = X] o '(0- ( 63 ) 

(0 (0 



For the background Eq. ([11]) remains valid for the collective fluid quantities. An additional equation follows from 
the individual energy-momentum conservation in Eq. ( ]6l| ) as 

+3if(/i(i) +P(i}) = Q(i)- (64) 
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For the perturbations Eqs. ( pL 2| - |1 8| ) remain valid for the collective fluid quantities. The additional equations we need 
follow from the individual energy-momentum conservation in Eq. (|6l|). From T^. b — Q^o and using Eq. (|l2] ) we 



have the energy conservation of the fluid components 

k 

5fi(i) + 3H + 8p {{) ) = (//(j) + A(i)« + (M(i) + P(i)) K + <^Q(i)- 



From k b — Q(i) Q we have the momentum conservation of the fluid components 



ff(l-&? ) 



1 + C 



(i) 



M(i) +P(i) 



fc 




W(t) = - 


a H f 


v a 


M(») +P(i) V 



2 /c 2 - 3K 



-<T(i) ~ J(i) 



(65) 



(66) 



By adding Eqs. (^,^6|) over the components we have Eqs. ( fi7|Jl8|) . In the multi-component situation Eqs. ( [l2f - 
|l8| , |65| , |66|) are the complete set expressed in a gauge- ready form. If we decompose the pressure as Sp^ = c^Sfi^ +e(i) 



with c 2 ^ = p\i)/fi(i), comparing with Eq. (|19|) we have 



e = 5p- c 2 Jn = ^2 [( c fi) ~ c l) <W) + e (i) 



Under the gauge transformation, from Eq. (20) in B, we have: 



(67) 



<*£(i) = <*M(i) ~ s P(i) = s P(i) - i>(i)t\ = v {i) £*, 



J, 



( - J M 



Thus, we have the following additional gauge conditions: 

<fyt (i ) = 0, 5p (i ) = 0, v {l )/k = 0, etc. 



(68) 
(69) 



Any one of these gauge conditions also fixes the temporal gauge condition completely. Studies of the multi-component 
situations can be found in [pU§J22|| . 



D. Rotation and Gravitational Wave 

For the vector-type perturbation, using Eqs. (p|JlC|), we have: 



SnGST^ = ± {a 2 + B J« + a (c^ + C*) ] }' , 

i(a^)o)- 



= -ST 



(70) 
(71) 
(72) 



where Eq. ( |72| ) follows from T*. 6 — 0, and apparently Eqs. (^,^l|) follow from the Einstein equations. We can 



show that Eq. (^) follows from Eqs. (|7(],|7l]). We introduce a gauge- invariant velocity related variable vi^x, t) as 

T^ v „ = (fi +p)Va . We can show Va <x auj where to is the amplitude of vorticity vector, see Sec. 5 of Thus, 
ignoring the anisotropic stress in the RHS of Eq. (|72) we have the conservation of angular momentum as 

Angular Momentum ~ a 3 (/x + p) x a x V^ w '(x, t) = constant in time, (73) 



which is valid for general K, A, and p(fi) in general scale. 

For the tensor-type perturbation, using Eqs. (^ 10) in the Einstein equations, we have 



(74) 



In the large-scale limit, assuming K = and ignoring the anisotropic pressure Eq. (|7J) has a general integral form 
solution 
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C$ (x, t) = C$ (x) -Dg(x) J ±dt, (75) 

where CVg and -D^g are integration constants for relatively growing and decaying solutions, respectively. Thus, ignoring 
the transient term the amplitude of the gravitational wave is conserved in the super-horizon scale considering general 
evolution of the background equation of state. 

The conservation of angular momentum and the equation for the gravitational wave were first derived in [[l]. 
Evolutions in the multi-component situation were considered in Sec. 5 of M. 



V. DISCUSSIONS 

We would like to make comments on related works in several textbooks: Eq. (15.10.57) in ]l2|] , Eq. (10.118) in p3[ , 
Eq. (11.5.2) in (24), Eq. (8.52) in |p5j , and problem 6.10 in ^(| are in error. All these errors are essentially about the 
same point involved with a fallible algebraic mistake in the synchronous gauge. The correction in the case of [12] ] was 
made in fl27j] : in a medium with a nonvanishing pressure the equation for the density fluctuation in the synchronous 
gauge becomes third order because of the presence of a gauge mode in addition to the physical growing and decaying 
solutions which is true even in the large-scale limit. The truncated second order equation i n |]l2| , p5| picks up a gauge 
mode instead of the physical decaying solution in the synchronous gauge. The errors in ]25| , |26| ] are based on imposing 
the synchronous gauge and the comoving gauge simultaneously, and thus happen to end up with the same truncated 
second order equation as in |l2| ] . In a medium with nonvanishing pressure one cannot impose the two gauge conditions 
simultaneously (even in the large-scale limit). In Sec. 11.5 of [£4| the authors proposed a simple modification of the 
Newtonian theory which again happens to end up with the same incorrect equation as in the other books. In the 
App endix we elaborate our point. For comments on other related errors in the literature, see Sec. 3.10 in p5[ and 



We would like to remark that these errors found in the synchronous gauge are not due to any esoteric aspect of 
the gauge choice in the relativistic perturbation analyses. Although not fundamental, the number of errors found in 
the literature concerning the synchronous gauge seem to indicate the importance of using the proper gauge condition 
in handling the problems. For our own choice of the preferred gauge-invariant variables suitable for handling the 



hydrodynamic perturbation and consequent analyses, see Sec. IV; the reasons for such choices are made in Sec. III. 
We wish to recall, however, that although one may need to do more algebra in tracing the remnant gauge mode, even 
the synchronous gauge condition is adequate for handling many problems as was carefully done in the original study 
by Lifshitz in 1946 fjjj. 

In this paper we have tried to identify the variables in the relativistic perturbation analysis which reproduce the 
correct Newtonian behavior in the prcssurclcss limit. In the first part we have shown that S in the comoving gauge (S v ) 
and v and ip in the zero-shear gauge (v x and ip x ) show the same behavior as the corresponding Newtonian variables in 
general scales. In fact, these results have already been presented in ]To| . Also, various general and as ym ptotic solutions 



wo 



for every variable in the pressureless medium are presented in the Tables of []1 Of - Compared with [|10fl , in Sec. Ill 
tried to reinforce the correspondence by explicitly showing the second order differential equations in several gauge 
conditions. We have also added some additional insights gained after publishing jjS. The second part contains some 
original results. Using the gauge- invariant variables S v , v x and ip x we write the relativistic hydrodynamic cosmological 
perturbation equations in Eqs. (f4l|-fi"4"|) ; actually, these equations are also known in [^J. The new results in this paper 
are the compact way of deriving the general large-scale solutions in Eqs. (]47j-[49"|) , and the clarification of the general 
large-scale conservation property of ip v in Eq. ( J53| ) . 

The underlying mathematical or physical reasons for the variables S v , ip x , v x , and tp v having the distinguished 
behaviors, compared with many other available gauge-invariant combinations, may still deserve further investigation. 
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APPENDIX: COMMON ERRORS IN THE SYNCHRONOUS GAUGE 



In the following we correct a minor confusion in the literature concerning perturbation analyses in the synchronous 
gauge. The argument is based on p7| . For simplicity we consider a situation with K = = A, e = = a and 
w = constant (thus c 2 = w). 

The equation for the density perturbation in the comoving gauge is given in Eq. (|45| ) . In our case we have 



S v + (2 - 3w)H6 v + 



k 2 

c 2 —- 4ttG^(1 - w)(l + 3w) 



0. 



(Al) 



The solution in the super-sound-horizon scale is presented in Eq. (|5C_ 

In the synchronous gauge, from Eqs. (|l6 - 18 ) we have (the subindex a indicates the synchronous gauge variable) 

,k 2 



5 a + 2HS a + 



c-— -47rG*Ai(l + w)(l + 3w) 



5 a + 3w(l + w) — Hv a = 0, 
a 



(A2) 



Va + (1 — 3w)Hv a 



a 1 



-5 a = 0. 



(A3) 



Notice that for w/0 these two equations are generally coupled even in the large-scale limit. From these two equations 
we can derive a third order differential equation for 5 a as 



11 -3w 



HS a 



24w - 9w z 



H 2 



,k 2 



4 V 



/) (1 + 3w) (-1 + 9w) H 3 + - (1 + w) He 2 — 



S a =0. 



In the large-scale limit the solutions are 



t 3 < 1 + "> , *3<l + w) , f 



(A4) 



(A5) 



In Appendix B of p7[ | we showed that the third solution with S a oc t 1 is nothing but a gauge mode for a medium 
with w ^ 0. [As mentioned before, the combination S a is not gauge-invariant. From Eq. ( pp[ ) we have 5 a = 
5 + 3H (1 + w) J* adt and the lower bound of the integration gives rise to the gauge mode which is proportional to 
t — 1 .] For a pressureless case the physical decaying solution also behaves as t^ 1 and the second solution in Eq. (A5) is 
invalid, see Sec. 4 in ^7j. In Eq. (B5) of [^7j we derived the relation between solutions in the two gauges explicitly; the 

2(l + 3w) 

growing solutions are the same in both gauges whereas the decaying solutions differ by a factor (k/aH) 2 oc 1 3 (!+ w ) . 

Now, we would like to point out that, in a medium with w ^ 0, one cannot ignore the last term in Eq. (A2) even in 
the large-scale limit. If we incorrectly ignore the last term in Eq. (A2) we recover the wrong equation in the textbooks 
mentioned in Sec. which is 

,k 2 



2H5 a + 



4ttG/j(1 + w)(1 + 3w) 



0. 



In the large-scale we have solutions 



2(l + 3w) 
t 3(l+w) 



(A6) 



(A7) 



By ignoring the last term in Eq. (A2), in the large-scale limit we happen to recover the fictitious gauge mode under 
the price of losing the physical decaying solution [the second solution in Eq. (A5)]. Thus, in the large-scale limit one 
cannot ignore the last term in Eq. (A2) in a medium with a general pressure; the reason is obvious if we see Eq. (A4). 
Also, one cannot impose both the synchronous gauge condition and the comoving gauge condition simultaneously. If 
we simultaneously impose such two conditions, thus setting v = = a, from Eq. ( |18| ) we have 







a 1 



A* 



2 k 2 



3K a 
— J, 



Thus, for e — = a and medium with the nonvanishing pressure we have <5 
argument remains valid even in the large-scale limit. 



(A8) 

which is a meaningless system; this 
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Table 1. Newtonian correspondences: For the synchronous gauge we ignore the gauge mode. Thus the syn- 
chronous gauge is equivalent to the comoving gauge. X indicates that the behavior differs from the Newtonian one. 
The small-scale implies the scale smaller than the (visual) horizon. Explicit forms of exact and asymptotic solutions 
for every variable are presented in [ ^o|JT^1 . 



Gauge 



Variable 



General Scale 



Small Scale 



Comoving gauge 


d v /k 


Newtonian 


Newtonian 


Zero-shear gauge 


S x 


X 


Newtonian 


Uniform-expansion gauge 




X 


Newtonian 


Synchronous gauge 


Sa 


Newtonian 


Newtonian 


Uniform-curvature gauge 


Sip 


X 


Newtonian 


Uniform-density gauge 


5 = 








Comoving gauge 


v = 








Zero-shear gauge 


v x 


Newtonian 


Newtonian 


Uniform-expansion gauge 


v K 


X 


Newtonian 


Synchronous gauge 


V a 








Uniform-curvature gauge 


v v 


X 


X 


Uniform-density gauge 


vs 


X 


X 


Comoving gauge 




X 


X 


Zero-shear gauge 


<Px 


Newtonian 


Newtonian 


Uniform-expansion gauge 




X 


Newtonian 


Synchronous gauge 




X 


X 


Uniform-curvature gauge 


Lf = 








Uniform-density gauge 




X 


X 



Note added: 

We introduce a variable 



K/a 2 



4ttG(h + p) 

From the general large-scale solutions in eqs. (47), (53) we can show 

$(x,i)=C(x), 



(B.l) 



(B.2) 



where the dominating decaying solution vanished. Thus, remarkably, $ is conserved, in the limit of vanishing c 2 s k 2 / a? 
term, considering general K, A, and time- varying p(/i). In a pressureless medium eq. (B.2) is an exact solution valid 
in general scale. 



1G 



